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We investigate the steady states and dynamical instabilities resulting from "particles" depositing 
on (fusion) and pinching off (fission) a fluid membrane. These particles could be either small lipid 
vesicles or isolated proteins. In the stable case, such fusion/fission events suppress long wavelength 
fluctuations of the membrane. In the unstable case, the membrane shoots out long tubular structures 
reminiscent of endosomal compartments or folded structures as in internal membranes like the 
endoplasmic reticulum or Golgi. 



PACS numbers: 87.f6.-b, 05.40.-a, 82.65.Dp 

The plasma membrane and internal membranes of eu- 
karyotic cells such as endosomal compartments, Golgi 
and endoplasmic reticulum (ER) are subject to random 
deposition and evaporation of 'particles' [Q. These par- 
ticles could be isolated lipids or membrane proteins but 
most often are small vesicles carrying cargo. For in- 
stance, in the endocytic pathway, small vesicles (diam- 
eter « 50nm), pinch off from the plasma membrane and 
fuse with the sorting endosome. The sorting endosome is 
seen to possess long, thin, tubular protrusions that spo- 
radically break away from it B. Likewise the ER < — > 
Golgi < — > Plasma Membrane system Q witnesses a ma- 
jor trafficking of protein carrying vesicles. The Golgi and 
the ER are highly ramified membranal structures, often 
possessing long thin tubular profusions |£|. 

The fusion/fission events are 'active' processes in that 
they are triggered by signalling proteins which recruit an 
elaborate array of proteins as part of the fusion/fission 
machinery pipl. These proteins, such as coat proteins, 
SNARES and rabs are functionally activated by energy 
provided by the consumption of ATP/GTP. Though at 
first sight the phenomenon appears alarmingly complex, 
we argue below that it gives rise to precise questions of a 
generic nature that do not require a knowledge of specific 
biomolecular details, and which may be addressed within 
the framework of nonequilibrium statistical physics. 

In this Letter we ask : Starting with an almost planar 
fluid membrane, what is its eventual shape when sub- 
ject to such active fusion/fission events ? We will find 
that the final fate of the membrane is (a) a steady state 
shape whose long wavelength fluctuations are suppressed 
by such events or (b) a dynamical instability giving rise to 
either long tubular structures or folded structures which 
grow rapidly. This rapid growth would be halted by other 
processes (not included in our analysis) which might lead 
to an eventual pinch off from the membrane. 

We take a coarse-grained approach in which the equi- 
librium shapes of a fluid membrane (in the absence of 
fusion/fission), parametrized by its 3-dimensional con- 
formation R(ui,U2) as a function of local coordinates 
(wi,M2), are governed by the local curvature 0] and the 



areal density -0 of 'particles' on it, through the effective 
free-energy functional 



Pel ^ ^ I (fu^ 



a + K,H + 2coV'-ff + XW 



(1) 



The local mean curvature H — 1/2 gab^ ■ dadb^, 
where N is the local unit (outward) normal and gab is 
the induced metric {g =det gab)- The density field tp 
has the following interpretation — when the particles 
are tiny lipid vesicles, ijj = p — po is the excess areal lipid 
density (/5 is the local areal lipid density and po is its 
equilbrium value) and cq — 0; when the particles are 
isolated membrane proteins or structurally asymmetric 
transmembrane proteins, ip is the areal protein density 
and the spontaneous curvature modulus cq ^ 0, reflect- 
ing the change in the local elastic properties of the mem- 
brane due to the embedded protein. 

In this paper we will work entirely in the Monge rep- 
resentation where the membrane shape is represented by 
a height field h{x, y) and all physical quantities are com- 
puted in the horizontal x = (a:, y) plane. In this repre- 
sentation g = 1 + (V/i)2 and H = \/^h/g^/^. The areal 
densities projected onto the x plane are obtained from 
P = P^/9- I^^ 0^1'' notation quantities with and with- 
out the tilde will denote intrinsic and projected densi- 
ties respectively. In this Monge representation, thermal 
equilibrium height correlations are given by (|ft.gp) c>c 
T/ (crg2 -I- KRq^), where kr = Kc- (kcCq) /x- 

We now start the fusion/fission processes which drive 
the membrane out of equilibrium, giving rise to a dy- 
namical excess density field V'(x: t) a-nd a dynamical 
shape /i(x, t). For convenience we shall use the pic- 
ture of fusion/fission of small vesicles in the rest of the 
paper, although our mathematical treatment includes 
the deposition and 'evaporation' of isolated proteins as 
well. Our subsequent analysis is valid over length scales 
much larger than the typical size £ of the 'particles' and 
time scales much longer than the duration of individual 
fusion(fission)-events r. Thus we shall ignore the compli- 
cated microscopic processes at work during the event. As 



in g , we try to extract the most robust features through 
general symmetry arguments. To facihtate analysis, we 
can regard a local region of the membrane as being in 
one of three states: (i) an undistorted membrane; (ii) a 
membrane with a 'bump'; (iii) an undistorted membrane 
with a proximate vesicle. The balance between (i) and 
(ii) is governed largely by elastic and equilibrium ther- 
modynamic forces whereas that between (ii) and (iii) is 
determined mainly by active processes M (Fig. 1). 
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FIG. 1. The "active" and "elastic" stages of (a) fusion 
and (b) fission. 

Define an 'fusion-event' e+(x, i) and an 'fission-event' 
e_(x, t) as random numbers taking values or 1, with a 
nonzero mean (e±) = e± and distributed identically and 
independently at different x and t. Every event gives rise 
to a change in the local height and density ; we would 
therefore like to calculate the rate of such fusion/fission 
events (number of events/area/time). We argue as fol- 
lows : Since the fusion/fission-events happen only on the 
cytosolic side in the case of internal membranes, and in 
an in-out asymmetric manner in the case of the plasma 
membrane, a local observer on the membrane can tell up 
from down in the Monge description. Thus all dynamical 
effects should reflect this up/down membrane asymme- 
try. This implies that the rate of fusion (fission) events, 
r"'"(r~), should be sensitive to the local mean curvature 
and hence to lowest order can be written as 



r+ = Ai + X2H + Xsp + XipH 
■Tn = -"iP + l^ipH + ... 



(2) 
(3) 



where the coefficients Ai, ^1 > by virtue of the defini- 
tion of the fusion/fission processes. Note that symmetry 
arguments say nothing about the values of the other pa- 
rameters, whose magnitudes and signs are determined by 
the microphysics of fusion/fission. 

To determine the net normal velocity, we multiply 
the above rates by a typical velocity (./t and include a 
permeative contribution proportional to the force f^j = 
<5Fei/<5R, 



{(./t) [e+(a;,t)r+ -e_(x,t)r; 



r/e. 



(4) 



where F is a mobility. Similarly the changes in the density 
due to the fusion/fission events are proportional to the 
"particle currents" , 



3+^{A+/t) e+(x,i)r+ 
j_ = (A_/r) e-(x,t)r- 



(5) 
(6) 



where A± are dimensionless, positive constants. It will 
turn out that we may consistently take A+ = A_ (for 



lipid vesicles, A^ 



/a, the amount of lipid material 



brought during a typical fusion event, a being the area 
per lipid). 

In the Monge representation, we decompose the normal 
velocity along the vertical z axis, v^ — Vni^ ■ z), and the 
horizontal x plane, w" = w„(N • x"), where a = 1,2 
labels the basis vectors in the x-y plane 0. Using the 
definitions (2) to (6) this yields the equations of motion, 

dth = - [e+ (Ai + AaV^/i + A3P + X^p^l'^h) 
-e_ [pip + p2P^'^h)\ 
-F [-crV^/i + KcV/i - KcCqVV] +Vh + fh (7) 

for the height field and 



dp 
'di 



-V.(pv^) + V.(^i5V^ 



+ j+ - i- + 



v-f, 



p 1 



(8) 



for 



/ 



the density 



field 



i. In Eq. (0), Vh = 

i2t^\2 ^ - TTlThrit)^ ^^ ^^^ vertical component of the 
viscous flow [rj is the solvent viscosity) induced at the 
membrane surface by the membrane elastic stresses [R. 
We have ignored other sources of dissipation coming from 
intramembrane viscosity and bilayer friction 0. The 
equations of motion also contain thermal noise sources 
having zero mean ; in the Monge representation, the vari- 
ance (|//i(q,<)P) = 2T T + {AT^\q\y^ , while the vari- 
ance of the vector noise fp is 2TD. Note the source and 
sink terms j± in Eq. (ph do not conserve particle material. 
Let us first check if Eqs. (|3,@) admit steady state so- 
lutions when we replace the noise terms by their means. 
We shall drop the hydrodynamic force for simplicity. Our 
subsequent analysis will be for a tensionless {a = 0) mem- 
brane with no spontaneous curvature (co — 0). We find 
spatially homogeneous steady states for which j-^ — j_, 
which implies that in the steady state the density attains 
a uniform value given by 



e+Ai 



-A, 



e-pi 



(9) 



while the steady state height is a constant (which we take 
to be 0). Since Ai > 0, the above steady state holds when 
B = e+Aa — e_^i < 0. Is this steady state linearly stable 
to perturbations in h and p = poo -I- </) ? In momentum 
space, the linearised perturbations are given by. 
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(10) 



where A = e+A2 + Poc{e+M — e_/i2)- The eigenvalues of 
the dynamical matrix for small q are 



Ai,: 



f,B-q^[Dx+iA] 

^cT-j^ADx/B 



ai) 



The nonconservative (j) ensures that one of the 
modes is always 'fast' i.e., has a nonzero relaxation 
rate at sufficiently long wavelength. If IA/t > 
-inm[Dx, Kcr{e'^/aT)\B\/Dx] = ~A*, then both 
modes decay and the homogeneous flat steady state is 
stable. We calculate the equal-time height fluctuations 
about this steady state, by first restoring the thermal and 
non-thermal noises £± — e± + de±, in Eq. (^^, fourier 
transforming with respect to x and t, and then integrat- 
ing the correlator with respect to frequency to to obtain, 



{h,it)h^g{t)) 



E 
7^ 



T 



M,2 



Kcr 



(12) 



for small q. The first term is a purely nonequilibrium 
contribution to the height fluctuations with a q^ disper- 
sion and an amplitude E related to the variance of Se± 
and the activity. The second term is a modification of 
the usual thermal correlator — the fission/fusion activ- 
ity suppresses height fluctuations at small q through an 
activity- induced tension (see Q) proportional to A. On 
the other hand, density correlations are short ranged with 
a correlation length ^ — B~^^^. 

If A < —A*, then either one or both modes are unsta- 
ble. The unstable massless mode grows as q'^ for small q. 
Physically this instability is a consequence of the curva- 
ture dependent rate of fusion : an initial outward bump 
(V^ft. < 0) will promote a faster rate of fusion, adding on 
more h. 

Going beyond linear stability analysis in the case A < 
—A* requires a numerical simulation of Eqs. (0|g). It is 
convenient to work in dimensionless variables, where x — > 
x/y'Kc/x, t -^ t/{Kc/Dx), h -^ h/e and p -^ pjilja^). 
Using an Euler discretisation on a square grid of size 
iV^ = 100^, we convert the differential equations to a 
coupled map. The spatial and temporal increments have 
been chosen to satisfy At < F^^ (Ax)'* and Ai < (Aa;)^ 
to ensure numerical stability. For most of our analysis 
we have fixed Ax = 0.1 and At = 10^^. We have taken 
open boundary conditions and have used a large enough 
system to ensure that no finite size effects corrupt our 
interpretation. We report numerical results only for the 
linearly unstable case. 

We start with the homogeneous steady state, /i(x, 0) — 
and p = Poo, and set the non-thermal noises to their 
average values, e± = e±. Introduce a small bump, a 
paraboloid of height Lq and width wq smoothly con- 
nected to the rest of the membrane. The peak height 



L[t) and the width w(t) grow initially, the width then 
saturates at q~^ ^ |A|~^/^, the length scale correspond- 
ing to the fastest growing mode determined from the 
linear analysis. The height now grows exponentially as 
L(t) ^ exp(A^t) (Fig. 2). Beyond this linear regime, 
whose time scale is set by A, the nonlinear term vp'S/^h 
(where v = e^X^^ — e_/i2) starts becoming comparable. 
The sign of v dictates the subsequent evolution ||l^. If 
;/ < 0, the nonlinear term accelerates the growth, produc- 
ing a long thin tubule (Fig. 2). Since the local curvature 
of the membrane becomes positive in the region where the 
main peak joins the rest of the membrane, this gives rise 
to side branches which dig into the membrane at a much 
slower rate (Fig. 2). The resulting morphology is a long 
tubule which grows very fast. This fast growth may even- 
tually be halted by molecular processes not contained in 
our coarse-grained model, leading to a fissioning of the 
tubule from the parent membrane. 
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FIG. 2. Tubule configuration (obtained by rotating the 
profile about the vertical axis) at t = 1.1 when v < Q. 
Inset shows the peak height L and the depth of the side 
branch D (dashed line is the initial exponential growth 
of L). The parameters v = —0.012 and A = —0.16. 

When ly > 0, the exponential rise of the peak height 
is cut off by a decrease in the rate j = j+ — j_ . At the 
same time the density at the peak diffuses out into the 
flanks in such a way as to compensate for the increase 
coming from the addition of 'particles'. This leads to a 
saturation of the density at the peak and asymptotically 
approaches ps = —ej^Xijv. The depth of the flanks in- 
creases exponentially fast till it is comparable with the 
slowly increasing peak height L. The resulting structure 
is a folded configuration as in Fig. 3. 

So far we had set the non-thermal noises e± to their 
average values ; we had therefore needed an initial bump 
on the membrane to seed the growth of tubules. The 
situation is unchanged if instead we let e± (x, f ) — 
e± -|- (5e±(x, t), where the fiuctuations are uncorrelated 



in space and time. In practice, this initial localised per- 
turbation about the steady state may be achieved by 
either phase separation of different lipid species consti- 
tuting the membrane [^ or noise. If we however set 
e±(x, i) = e±{t)S{x), where e±(i) is stochastic, then the 
membrane spontaneously creates a bump at the fusion 
site, growing exponentially fast. 
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FIG. 3. Folded configuration (obtained on rotation 
about the vertical) at t = 2.6 when v > Q. Inset shows 
the height of the central peak L and the largest depth of 
the folds D (dashed line is the initial exponential growth 
of L). The parameters v = 0.012 and A = -0.14. 

The tubular extensions produced when // < look sim- 
ilar to the structures exhibited by the sorting endosome 
associated with the endocytic pathway of transferrin and 
LDL proteins. The sorting endosome is seen to possess 
several tubular structures of width « 50nm which grow 
to a length of about a micron before peeling off from 
the parent structure ||]. On the other hand, the folded 
structures produced when v > are most dramatically 
exhibited in a closed membrane, which may be studied 
using a dynamical triangulation monte carlo modified to 
add and subtract particles ||ll| . 

Such tubulation or folded morphologies could also be 
generated if the 'particles' depositing on and 'evaporat- 
ing' from the membrane are isolated proteins. The rates 
of these events may be altered by the action of some 
drug, such as Brefeldin-A. Our analysis suggests that 
such deposition-evaporation of proteins from a membrane 
could result in tubulation ll^ ]. 

A second look at Eqs. ([jj^) ^Uows us to make some def- 
inite predictions. We see that the surface tension, which 
is the coefficient of the V^ft. term, gets renormalized by 
an 'activity induced' factor, a = a -\- £ A/T. This renor- 
malization is also reflected in the form of the correlators, 
Eq. ( |l2| ) . Alternatively, applying a tension a to the mem- 
brane will affect the rate of fusion/fission events ll3]. 

The activity induced tension is in general inhomoge- 
ncous since the 'activity' A may be space dependent. 
This gives rise to gradients in the activity induced sur- 
face tension, driving flow of membrane constituents to- 



wards regions of lower tension. Ignoring bilayer friction, 
the velocity of this Marangoni flow may be written as 

vj^ — — — V(7, where ^ is the in-plane shear viscosity and 
^ is a length scale associated with the 'particle' size. The 
direction of flow is such that the membrane constituents 
will move towards the tubular extensions produced as a 
result of the activity. 

In future we would like to study the effects of stochastic 
fusion/fission of 'particles' on a closed fluid membrane 
|ll| , and to explore the complete nonequilibrium steady 
state diagram by solving the stochastic equations (n,H). 
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